Integrating taxonomic signals from MAGs and contigs improves read annotation and taxonomic profiling of metagenomes.

Metagenomic analysis typically includes read-based taxonomic profiling, assembly, and binning of metagenome-assembled genomes (MAGs). Here we integrate these steps in Read Annotation Tool (RAT), which uses robust taxonomic signals from MAGs and contigs to enhance read annotation. RAT reconstructs taxonomic profiles with high precision and sensitivity, outperforming other state-of-the-art tools. In high-diversity groundwater samples, RAT annotates a large fraction of the metagenomic reads, calling novel taxa at the appropriate, sometimes high taxonomic ranks. Thus, RAT integrative profiling provides an accurate and comprehensive view of the microbiome from shotgun metagenomics data. The package of Contig Annotation Tool (CAT), Bin Annotation Tool (BAT), and RAT is available at https://github.com/MGXlab/CAT_pack (from CAT pack v6.0). The CAT pack now also supports Genome Taxonomy Database (GTDB) annotations.

A low-cost open-source imaging platform reveals spatiotemporal insight into leaf elongation and movement.

Plant organs move throughout the diurnal cycle, changing leaf and petiole positions to balance light capture, leaf temperature, and water loss under dynamic environmental conditions. Upward movement of the petiole, called hyponasty, is one of several traits of the shade avoidance syndrome (SAS). SAS traits are elicited upon perception of vegetation shade signals such as far-red light (FR) and improve light capture in dense vegetation. Monitoring plant movement at a high temporal resolution allows studying functionality and molecular regulation of hyponasty. However, high temporal resolution imaging solutions are often very expensive, making this unavailable to many researchers. Here, we present a modular and low-cost imaging set-up, based on small Raspberry Pi computers, that can track leaf movements and elongation growth with high temporal resolution. We also developed an open-source, semi-automated image analysis pipeline. Using this setup we followed responses to FR enrichment, light intensity, and their interactions. Tracking both elongation and the angle of the petiole, lamina, and entire leaf in Arabidopsis (Arabidopsis thaliana) revealed insight into R:FR sensitivities of leaf growth and movement dynamics and the interactions of R:FR with background light intensity. The detailed imaging options of this system allowed us to identify spatially separate bending points for petiole and lamina positioning of the leaf.

Radiation protection in CT-guided interventions: does real-time dose visualisation lead to a reduction in radiation dose to participating radiologists? A single-centre evaluation.

AIM: To evaluate if real-time dose visualisation during computed tomography (CT)-guided interventions leads to a reduction in radiation dose to participating radiologists. MATERIALS AND METHODS: The individual radiation dose radiologists are exposed to during CT interventions was measured using dedicated dosimeters (RaySafe i2-system, Unfors RaySafe GmbH, Billdal, Sweden) worn over the usual radiation protective apron. Initially, only the total radiation dose was measured, without visualisation (control group). In the following study period, the radiation dose was shown to participants on a live screen in real-time (experimental group). In both groups, the dose was recorded in 1-second intervals. The results collected were evaluated by comparison using descriptive statistics and mixed-effect models. In particular, the variables experience, gender, role, and position during the intervention were analysed. RESULTS: In total, 517 measurements of 304 interventions (n=249 with and n=268 without live screen) performed by 29 radiologists acting as interventionalists or assistants were analysed. All CT-guided interventions were performed percutaneously, the majority of which (n=280) were microwave ablations (MWA). Radiation doses in the group without visualisation were comparable with usual dose rates for the corresponding intervention type. The mean total radiation dose was reduced by 58.1% (11.6 versus 4.86 µSv) in the experimental group (p=0.034). The highest reduction of 78.5% (15.55 versus 3.35 µSv) was observed in radiologists with the role of assistant (p=0.002). Sub-analysis showed significant dose reduction (p<0.0001) for the use of live screen in general; considering all variables, the role "assistant" alone had a statistically significant influence (p=0.002). CONCLUSION: The real-time visualisation of active radiation dose during CT interventions leads to a relevant reduction in radiation dose to participating radiologists.

The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population.

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

Cryptic genetic variation of expression quantitative trait locus architecture revealed by genetic perturbation in Caenorhabditis elegans.

Genetic perturbation in different genetic backgrounds can cause a range of phenotypes within a species. These phenotypic differences can be the result of the interaction between the genetic background and the perturbation. Previously, we reported that perturbation of gld-1, an important player in the developmental control of Caenorhabditis elegans, released cryptic genetic variation (CGV) affecting fitness in different genetic backgrounds. Here, we investigated the change in transcriptional architecture. We found 414 genes with a cis-expression quantitative trait locus (eQTL) and 991 genes with a trans-eQTL that were specifically found in the gld-1 RNAi treatment. In total, we detected 16 eQTL hotspots, of which 7 were only found in the gld-1 RNAi treatment. Enrichment analysis of those 7 hotspots showed that the regulated genes were associated with neurons and the pharynx. Furthermore, we found evidence of accelerated transcriptional aging in the gld-1 RNAi-treated nematodes. Overall, our results illustrate that studying CGV leads to the discovery of hidden polymorphic regulators.

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds.

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Local light signaling at the leaf tip drives remote differential petiole growth through auxin-gibberellin dynamics.

Although plants are immobile, many of their organs are flexible to move in response to environmental cues. In dense vegetation, plants detect neighbors through far-red light perception with their leaf tip. They respond remotely, with asymmetrical growth between the abaxial and adaxial sides of the leafstalk, the petiole. This results in upward movement that brings the leaf blades into better lit zones of the canopy. The plant hormone auxin is required for this response, but it is not understood how non-differential leaf tip-derived auxin can remotely regulate movement. Here, we show that remote signaling of far-red light promotes auxin accumulation in the abaxial petiole. This local auxin accumulation is facilitated by reinforcing an intrinsic directionality of the auxin transport protein PIN3 on the petiole endodermis, as visualized with a PIN3-GFP line. Using an auxin biosensor, we show that auxin accumulates in all cell layers from endodermis to epidermis in the abaxial petiole, upon far-red light signaling in the remote leaf tip. In the petiole, auxin elicits a response to both auxin itself as well as a second growth promoter; gibberellin. We show that this dual regulation is necessary for hyponastic leaf movement in response to light. Our data indicate that gibberellin is required to permit cell growth, whereas differential auxin accumulation determines which cells can grow. Our results reveal how plants can spatially relay information about neighbor proximity from their sensory leaf tips to the petiole base, thus driving adaptive growth.

Coumarin biosynthesis genes are required after foliar pathogen infection for the creation of a microbial soil-borne legacy that primes plants for SA-dependent defenses.

Plants deposit photosynthetically-fixed carbon in the rhizosphere, the thin soil layer directly around the root, thereby creating a hospitable environment for microbes. To manage the inhabitants of this nutrient-rich environment, plant roots exude and dynamically adjust microbe-attracting and -repelling compounds to stimulate specific members of the microbiome. Previously, we demonstrated that foliar infection of Arabidopsis thaliana by the biotrophic downy mildew pathogen Hyaloperonospora arabidopsidis (Hpa) leads to a disease-induced modification of the rhizosphere microbiome. Soil conditioned with Hpa-infected plants provided enhanced protection against foliar downy mildew infection in a subsequent population of plants, a phenomenon dubbed the soil-borne legacy (SBL). Here, we show that for the creation of the SBL, plant-produced coumarins play a prominent role as coumarin-deficient myb72 and f6'h1 mutants were defective in creating a Hpa-induced SBL. Root exudation profiles changed significantly in Col-0 upon foliar Hpa infection, and this was accompanied by a compositional shift in the root microbiome that was significantly different from microbial shifts occurring on roots of Hpa-infected coumarin-deficient mutants. Our data further show that the Hpa-induced SBL primes Col-0 plants growing in SBL-conditioned soil for salicylic acid (SA)-dependent defenses. The SA-signaling mutants sid2 and npr1 were unresponsive to the Hpa-induced SBL, suggesting that the protective effect of the Hpa-induced shift in the root microbiome results from an induced systemic resistance that requires SA-signaling in the plant.

A metagenomic portrait of the microbial community responsible for two decades of bioremediation of poly-contaminated groundwater.

Biodegradation of pollutants is a sustainable and cost-effective solution to groundwater pollution. Here, we investigate microbial populations involved in biodegradation of poly-contaminants in a pipeline for heavily contaminated groundwater. Groundwater moves from a polluted park to a treatment plant, where an aerated bioreactor effectively removes the contaminants. While the biomass does not settle in the reactor, sediment is collected afterwards and used to seed the new polluted groundwater via a backwash cycle. The pipeline has successfully operated since 1999, but the biological components in the reactor and the contaminated park groundwater have never been described. We sampled seven points along the pipeline, representing the entire remediation process, and characterized the changing microbial communities using genome-resolved metagenomic analysis. We assembled 297 medium- and high-quality metagenome-assembled genome sequences representing on average 46.3% of the total DNA per sample. We found that the communities cluster into two distinct groups, separating the anaerobic communities in the park groundwater from the aerobic communities inside the plant. In the park, the community is dominated by members of the genus Sulfuricurvum, while the plant is dominated by generalists from the order Burkholderiales. Known aromatic compound biodegradation pathways are four times more abundant in the plant-side communities compared to the park-side. Our findings provide a genome-resolved portrait of the microbial community in a highly effective groundwater treatment system that has treated groundwater with a complex contamination profile for two decades.

Disentangling the genetic basis of rhizosphere microbiome assembly in tomato.

Microbiomes play a pivotal role in plant growth and health, but the genetic factors involved in microbiome assembly remain largely elusive. Here, we map the molecular features of the rhizosphere microbiome as quantitative traits of a diverse hybrid population of wild and domesticated tomato. Gene content analysis of prioritized tomato quantitative trait loci suggests a genetic basis for differential recruitment of various rhizobacterial lineages, including a Streptomyces-associated 6.31 Mbp region harboring tomato domestication sweeps and encoding, among others, the iron regulator FIT and the water channel aquaporin SlTIP2.3. Within metagenome-assembled genomes of root-associated Streptomyces and Cellvibrio, we identify bacterial genes involved in metabolism of plant polysaccharides, iron, sulfur, trehalose, and vitamins, whose genetic variation associates with specific tomato QTLs. By integrating 'microbiomics' and quantitative plant genetics, we pinpoint putative plant and reciprocal rhizobacterial traits underlying microbiome assembly, thereby providing a first step towards plant-microbiome breeding programs.

The genetic architecture underlying body-size traits plasticity over different temperatures and developmental stages in Caenorhabditis elegans.

Most ectotherms obey the temperature-size rule, meaning they grow larger in a colder environment. This raises the question of how the interplay between genes and temperature affects the body size of ectotherms. Despite the growing body of literature on the physiological life-history and molecular genetic mechanism underlying the temperature-size rule, the overall genetic architecture orchestrating this complex phenotype is not yet fully understood. One approach to identify genetic regulators of complex phenotypes is quantitative trait locus (QTL) mapping. Here, we explore the genetic architecture of body-size phenotypes, and plasticity of body-size phenotypes at different temperatures using Caenorhabditis elegans as a model ectotherm. We used 40 recombinant inbred lines (RILs) derived from N2 and CB4856, which were reared at four different temperatures (16, 20, 24, and 26 °C) and measured at two developmental stages (L4 and adult). The animals were measured for body length, width at vulva, body volume, length/width ratio, and seven other body-size traits. The genetically diverse RILs varied in their body-size phenotypes with heritabilities ranging from 0.0 to 0.99. We detected 18 QTL underlying the body-size traits across all treatment combinations, with the majority clustering on Chromosome X. We hypothesize that the Chromosome X QTL could result from a known pleiotropic regulator-npr-1-known to affect the body size of C. elegans through behavioral changes. We also found five plasticity QTL of body-size traits where three colocalized with body-size QTL. In conclusion, our findings shed more light on multiple loci affecting body-size plasticity and the possibility of co-regulation of traits and traits plasticity by the same loci under different environments.

2D morphometric analysis of Arabidopsis thaliana nuclei reveals characteristic profiles of different cell types and accessions.

Functional changes of cells upon developmental switches and in response to environmental cues are often reflected in nuclear phenotypes, showing distinctive chromatin states corresponding to transcriptional changes. Such characteristic nuclear shapes have been microscopically monitored and can be quantified after differential staining of euchromatin and heterochromatin domains. Here, we examined several nuclear parameters (size, DNA content, DNA density, chromatin compaction, relative heterochromatin fraction (RHF), and number of chromocenters) in relation to spatial distribution of genes and transposon elements (TEs), using standard 2D fluorescence microscopy. We provide nuclear profiles for different cell types and different accessions of Arabidopsis thaliana. A variable, yet significant, fraction of TEs was found outside chromocenters in all cell types, except for guard cells. The latter cell type features nuclei with the highest level of chromatin compaction, while their chromocenters seem to contain gene-rich regions. The highest number of parameter correlations was found in the accession Cvi, whereas Ler showed only few correlations. This may point at differences in phenotype robustness between accessions. The significantly high association of NOR chromocenters in accessions Ws and Cvi corresponds to their low RHF level.

Embodiment in distributed information processing: "Solid" plants versus "liquid" ant colonies.

Information processing is an essential part of biology, enabling coordination of intra-organismal processes such as development, environmental adaptation and inter-organismal communication. Whilst in animals with specialised brain tissue a substantial amount of information processing occurs in a centralised manner, most biological computing is distributed across multiple entities, such as cells in a tissue, roots in a root system or ants in a colony. Physical context, called embodiment, also affects the nature of biological computing. While plants and ant colonies both perform distributed computing, in plants the units occupy fixed positions while individual ants move around. This distinction, solid versus liquid brain computing, shapes the nature of computations. Here we compare information processing in plants and ant colonies, highlighting how similarities and differences originate in, as well as make use of, the differences in embodiment. We end with a discussion on how this embodiment perspective may inform the debate on plant cognition.

Mechanisms of far-red light-mediated dampening of defense against Botrytis cinerea in tomato leaves.

Plants detect neighboring competitors through a decrease in the ratio between red and far-red light (R:FR). This decreased R:FR is perceived by phytochrome photoreceptors and triggers shade avoidance responses such as shoot elongation and upward leaf movement (hyponasty). In addition to promoting elongation growth, low R:FR perception enhances plant susceptibility to pathogens: the growth-defense tradeoff. Although increased susceptibility in low R:FR has been studied for over a decade, the associated timing of molecular events is still unknown. Here, we studied the chronology of FR-induced susceptibility events in tomato (Solanum lycopersicum) plants pre-exposed to either white light (WL) or WL supplemented with FR light (WL+FR) prior to inoculation with the necrotrophic fungus Botrytis cinerea (B.c.). We monitored the leaf transcriptional changes over a 30-h time course upon infection and followed up with functional studies to identify mechanisms. We found that FR-induced susceptibility in tomato is linked to a general dampening of B.c.-responsive gene expression, and a delay in both pathogen recognition and jasmonic acid-mediated defense gene expression. In addition, we found that the supplemental FR-induced ethylene emissions affected plant immune responses under the WL+FR condition. This study improves our understanding of the growth-immunity tradeoff, while simultaneously providing leads to improve tomato resistance against pathogens in dense cropping systems.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population.

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5-46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype-phenotype relationship.

Transcriptome Signatures in Pseudomonas simiae WCS417 Shed Light on Role of Root-Secreted Coumarins in Arabidopsis-Mutualist Communication.

Pseudomonas simiae WCS417 is a root-colonizing bacterium with well-established plant-beneficial effects. Upon colonization of Arabidopsis roots, WCS417 evades local root immune responses while triggering an induced systemic resistance (ISR) in the leaves. The early onset of ISR in roots shows similarities with the iron deficiency response, as both responses are associated with the production and secretion of coumarins. Coumarins can mobilize iron from the soil environment and have a selective antimicrobial activity that impacts microbiome assembly in the rhizosphere. Being highly coumarin-tolerant, WCS417 induces the secretion of these phenolic compounds, likely to improve its own niche establishment, while providing growth and immunity benefits for the host in return. To investigate the possible signaling function of coumarins in the mutualistic Arabidopsis-WCS417 interaction, we analyzed the transcriptome of WCS417 growing in root exudates of coumarin-producing Arabidopsis Col-0 and the coumarin-biosynthesis mutant f6'h1. We found that coumarins in F6'H1-dependent root exudates significantly affected the expression of 439 bacterial genes (8% of the bacterial genome). Of those, genes with functions related to transport and metabolism of carbohydrates, amino acids, and nucleotides were induced, whereas genes with functions related to cell motility, the bacterial mobilome, and energy production and conversion were repressed. Strikingly, most genes related to flagellar biosynthesis were down-regulated by F6'H1-dependent root exudates and we found that application of selected coumarins reduces bacterial motility. These findings suggest that coumarins' function in the rhizosphere as semiochemicals in the communication between the roots and WCS417. Collectively, our results provide important novel leads for future functional analysis of molecular processes in the establishment of plant-mutualist interactions.

The chemical compound 'Heatin' stimulates hypocotyl elongation and interferes with the Arabidopsis NIT1-subfamily of nitrilases.

Temperature passively affects biological processes involved in plant growth. Therefore, it is challenging to study the dedicated temperature signalling pathways that orchestrate thermomorphogenesis, a suite of elongation growth-based adaptations that enhance leaf-cooling capacity. We screened a chemical library for compounds that restored hypocotyl elongation in the pif4-2-deficient mutant background at warm temperature conditions in Arabidopsis thaliana to identify modulators of thermomorphogenesis. The small aromatic compound 'Heatin', containing 1-iminomethyl-2-naphthol as a pharmacophore, was selected as an enhancer of elongation growth. We show that ARABIDOPSIS ALDEHYDE OXIDASES redundantly contribute to Heatin-mediated hypocotyl elongation. Following a chemical proteomics approach, the members of the NITRILASE1-subfamily of auxin biosynthesis enzymes were identified among the molecular targets of Heatin. Our data reveal that nitrilases are involved in promotion of hypocotyl elongation in response to high temperature and Heatin-mediated hypocotyl elongation requires the NITRILASE1-subfamily members, NIT1 and NIT2. Heatin inhibits NIT1-subfamily enzymatic activity in vitro and the application of Heatin accordingly results in the accumulation of NIT1-subfamily substrate indole-3-acetonitrile in vivo. However, levels of the NIT1-subfamily product, bioactive auxin (indole-3-acetic acid), were also significantly increased. It is likely that the stimulation of hypocotyl elongation by Heatin might be independent of its observed interaction with NITRILASE1-subfamily members. However, nitrilases may contribute to the Heatin response by stimulating indole-3-acetic acid biosynthesis in an indirect way. Heatin and its functional analogues present novel chemical entities for studying auxin biology.

Soil microbial legacies differ following drying-rewetting and freezing-thawing cycles.

Climate change alters frequencies and intensities of soil drying-rewetting and freezing-thawing cycles. These fluctuations affect soil water availability, a crucial driver of soil microbial activity. While these fluctuations are leaving imprints on soil microbiome structures, the question remains if the legacy of one type of weather fluctuation (e.g., drying-rewetting) affects the community response to the other (e.g., freezing-thawing). As both phenomenons give similar water availability fluctuations, we hypothesized that freezing-thawing and drying-rewetting cycles have similar effects on the soil microbiome. We tested this hypothesis by establishing targeted microcosm experiments. We created a legacy by exposing soil samples to a freezing-thawing or drying-rewetting cycle (phase 1), followed by an additional drying-rewetting or freezing-thawing cycle (phase 2). We measured soil respiration and analyzed soil microbiome structures. Across experiments, larger CO2 pulses and changes in microbiome structures were observed after rewetting than thawing. Drying-rewetting legacy affected the microbiome and CO2 emissions upon the following freezing-thawing cycle. Conversely, freezing-thawing legacy did not affect the microbial response to the drying-rewetting cycle. Our results suggest that drying-rewetting cycles have stronger effects on soil microbial communities and CO2 production than freezing-thawing cycles and that this pattern is mediated by sustained changes in soil microbiome structures.

Balancing Selection of the Intracellular Pathogen Response in Natural Caenorhabditis elegans Populations.

Genetic variation in host populations may lead to differential viral susceptibilities. Here, we investigate the role of natural genetic variation in the Intracellular Pathogen Response (IPR), an important antiviral pathway in the model organism Caenorhabditis elegans against Orsay virus (OrV). The IPR involves transcriptional activity of 80 genes including the pals-genes. We examine the genetic variation in the pals-family for traces of selection and explore the molecular and phenotypic effects of having distinct pals-gene alleles. Genetic analysis of 330 global C. elegans strains reveals that genetic diversity within the IPR-related pals-genes can be categorized in a few haplotypes worldwide. Importantly, two key IPR regulators, pals-22 and pals-25, are in a genomic region carrying signatures of balancing selection, suggesting that different evolutionary strategies exist in IPR regulation. We infected eleven C. elegans strains that represent three distinct pals-22 pals-25 haplotypes with Orsay virus to determine their susceptibility. For two of these strains, N2 and CB4856, the transcriptional response to infection was also measured. The results indicate that pals-22 pals-25 haplotype shapes the defense against OrV and host genetic variation can result in constitutive activation of IPR genes. Our work presents evidence for balancing genetic selection of immunity genes in C. elegans and provides a novel perspective on the functional diversity that can develop within a main antiviral response in natural host populations.